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We demonstrate the role of selectivity variation in the structure of the non-equilibrium extended 
space-charge using ID analytic and 2D numerical Poisson-Nernst-Planck models for the electro- 
diffusive transport of a symmetric electrolyte. This provides a deeper understanding of the underly¬ 
ing mechanism behind a previously-observed maximum in the resistance-voltage curve for a shallow 
micro-nanochannel interface device [Schiffbauer, Liel, Leibowitz, Park, and Yossifon, submitted to 
Phys. Rev. E.]. The current study helps to establish a connection between parameters such as 
the geometry and nanochannel surface-charge and the control of selectivity and resistance in the 
over-limiting current regime. 
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I. INTRODUCTION 

The current-voltage response of fabricated micro¬ 
nanochannel devices and ion-selective membranes are 
in many respects similar. Both exhibit ion-selective 
transport and produce concentration polarization (CP) 
in adjacent solution resulting in a near limiting sat¬ 
uration of the DC current response. At sufficiently 
high voltages, both exhibit an over-limiting current 
(OLC) which may be associated with electro-convective 
vortices among other mechanisms However, 

there are a number of differences between ion-selective 
membranes, such as Nafion, and fabricated nanochannel 
devices. The geometry of pores within a membrane are 
often irregularly shaped, highly tortuous, and poorly 
connected and the fixed charge density quite high. 
One critical difference is that, while membranes such as 
Nafion typically exhibit near ideal selectivity across a 
wide range of parameters [ToJ, nanochannel selectivity 
can vary appreciably. The similarity between membranes 
and nanochannels implies that a non-ideal membrane 
model can serve as a simple model for electro-diffusive 
transport through a nanochannel. 

Previous studies have shown variation in selectivity 
in the under-limiting regime El and demonstrated the 
role of non-ideality in rectification due to diffusion-layer 
asymmetry [l2|. Other studies concern the role of 
a non-equilibrium extended space charge (ESC) on 
transport of competing counter-ions through an ideal 
membrane fl3| . or consider space-charge dynamics for 
a binary electrolyte through a non-ideal membrane in 
the over-limiting current (OLC) regime ITij . However, 
these previous studies omit the first-order correction to 
the salt-concentration in the space charge regime and 
its coupling to the selectivity. Here we demonstrate the 
necessity of including these terms in capturing the effects 
of changing selectivity on the resistance maximum. This 
is expected to have consequences for the observability 
of the resistance maximum in real systems, where other 
OLC mechanisms can enhance or compete with the ESC 
driven resistance maximum [h}. Furthermore, it may 


have implications for the control of selectivity and ESC 
structure in practical applications, such as bio-molecular 
sensing, analyte pre-concentration, or electrokinetic 
desalination. 


II. THEORY 


In the following, we consider one-dimensional steady- 
state electro-diffusive transport of a symmetric, binary 
electrolyte through a system consisting of two solution 
layers flanking a non-ideally cation-selective membrane 
(see Fig. [TJ. We omit effects such as surface conduction 
and fluid-flow for simplicity. In particular, we focus our 
attention on the modulation of charge-selectivity in the 
over-limiting regime and its effect on the previously ob¬ 
served resistance maximum [n|. Using the steady-state 
(non-dimensional) Poisson-Nernst-Planck (PNP) equa¬ 
tions as a starting point, the ion concentration field is 
given by, 


V(x) 


dc± d(j) 

■*^ c±(l) s 


= J± 


(1) 


for each species ±, and the Poisson equation describes 
the electric potential 

- e 2 ^ = c+ (x) - c- ( x ) - N(x) (2) 

where the (negative) fixed charge density is N(x) = N 
for x £ {—1,1} and and 0 elsewhere. Lengths have been 
scaled by the membrane half-thickness, l, concentration 
by the bulk concentration, c D , fixed at the ends, and the 
electric potential by the thermal voltage, kT/e = 0.0254 
V. The ion diffusivity is scaled by the bulk value so that 
T>{x) = 1 in the electrolyte, but in general may have a 
different value, D m , in the membrane interior. The ap¬ 
plied field, hence electric current, is directed towards the 
right (positive) so the ESC sub-layer forms at the x = — 1 
interface as indicated in Fig. [I] The dimensionless Debye 
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FIG. 1: Basic geometry for the ID model. Regions a and d are 
respectively the depleted and enriched quasi-electroneutral 
diffusion layers, region c is charge-selective, and region b ap¬ 
proximates the width of the ESC. 


length, e = \Jee 0 kT / e 2 Lc 0 is taken as a small parameter, 
E<1. Its appearance in Eqn. [5] multiplying the highest 
derivative of the potential introduces the possibility of 
boundary layers and singular perturbation methods are 
typically used to obtain an analytic solution to the full 
problem (T^ |. 

The solutions to the ID problem are well-known 
in various approximations. Here we review results rel¬ 
evant to the non-ideal case. The problem is re-cast in 
terms of c = (c+ + c_)/2 (salt concentration,) the salt- 
flux density J+ = (j+ +j-)/2, and charge flux density 
J_ = (j + — j-)/2. The two flux densities are related by 
the selectivity factor G = J_ / J+ and the electric current 
density is I = 2 J_. While the full solution requires ei¬ 
ther singular perturbation theory or numerical methods, 
a useful first approximation may be obtained by taking 
the leading-order solution, or equivalently, setting e = 0, 
which is tantamount to imposing local electro-neutrality 
(LEN). These locally electro-neutral solutions for the salt 
concentration in the ID problem are bi-linear functions 
of the salt-flux density and the position, 

ci{x) = 1 — J+(L + x), (3) 

in the depleted (left) side and 

Cr{x) = J+{L - x) + 1 (4) 

in the enriched (right) side. The potentials are given by 
the integrals over the inverse concentration, or 

4>{x) = Gin [1 — J+(L + a;)] (5) 

for the depleted region and 

</>(x) = G\n[J+{L-x) + l]-V (6) 

for the potential in the enriched region. The difficulty 
with the LEN approximation can be seen clearly in 
Eqns. [3] and [5] As the current approaches a value cor¬ 
responding to the salt-flux density J+^um = 1 /(L — 1), 
the salt concentration at the depleted interface, x = — 1, 
approaches zero and the corresponding electric potential 
tends to infinity, i.e. it requires an infinite voltage drop 
to reach this value of current. However this current can 
not only be reached but exceeded in real systems. De¬ 
pending on the details of the system [§, 8j, the mecha¬ 
nisms which drive OLC can depend on violations of local 


electro-neutrality via the formation of the ESC [13- 
Inside the ion-selective region, x £ {—1,1}, we assume 
local quasi-electro-neutrality to hold to arbitrarily high 
currents. So Eqn. [2] reduces to c+ = c_ + N. The ques¬ 
tion of practical limits on this condition, especially in the 
context of the micro-nanochannel interface, is beyond the 
scope of the present study. The governing equations for 
the salt concentration and potential in the membrane in¬ 
terior are obtained from Eqns. [l] and [H 


and 


dc m N d(j) m _ J+ 
dx 2 dx V m 


dtfirn J— — 1 


( 7 ) 

( 8 ) 


These are integrated to yield the following transcendental 
equation for the concentration profile, 

[c m (x) - c m (-l)] + f Gin - 

and 

<!>m{x) =-^-${x) + 4> m {- 1) (10) 

I'm 

for the electric potential inside the ion selective region, 
where the integral 

SO) = J ^ c m( x ') dx ' (ID 


,(*)-- 




(9) 


is defined for convenience. 

For salt-flux densities lower than the limiting value, the 
response of the system may be obtained from the bound¬ 
ary conditions at x = ±L by assuming continuity of the 
electrochemical potentials across the quasi-equilibrium 
electric double layers (EDLs) between each region, as 
first proposed by Kirkwood [18|. Strictly speaking, this 
is an assumption of local, homogenous thermodynamic 
equilibrium between two physically distinct phases, sep¬ 
arated by a defining, infinitesimally thin Gibbs dividing 
plane. Because the EDLs are quite thin and the struc¬ 
ture appears to maintain a quasi-equilibrium character 
under current, this approximation may be employed in 
non-equilibrium scenarios. The extent of the validity of 
the internal quasi-equilibrium structure has been exam¬ 
ined thoroughly and shown to be maintained well into 
the non-equilibrium regime [16] . The system retains a 
quasi-equilibrium EDL as the innermost sub-layer, with 
the ESC developing as an adjacent structure [1J]. Thus 
one should always be able to impose, at least approxi¬ 
mately, a continuity condition to jump across this inner 
EDL boundary, resulting in a Donnan-like equilibrium 
between the interior of the ESC (adjacent to the thin 
EDL) and the membrane interior. We employ this as¬ 
sumption herein. 
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For the full solution, we follow the usual procedure 
for obtaining a master equation for the scaled E-field, 
E = —ed(j)/dx (see for instance Ref p,6|.) The main re¬ 
sults are the following equation for the scaled E-field, 

d? E E ^ 

e 2 ^- — + 2 J + E( x -x°) = -2eGJ + (12) 
and the equation for the cation concentration, 

c+{x) = ^E' + E - J + (x-x 0 ) (13) 

and the salt concentration, 

E% 

c{x) = — - J + ( x - x °) (14) 

Thus the concentrations may be calculated once a so¬ 
lution for the E-field is obtained. Note the integration 
constant x 0 . It can be associated with the intersection of 
the extrapolated depleted linear salt concentration profile 
to zero concentration. At and above the limiting current, 
this point moves into the region x € {—L, —1}, yielding 
an estimate of the edge position of the ESC. 

Asymptotic expansions for the scaled E-field are given 
by the following for regions sufficiently far (^> e 2//3 ) from 
either side of the point x = x Q : 


E{x) 


-eG 
( x—x 0 ) 


6 3 (4G-G 3 ) 
4J_|_(x— x 0 ) A 


0(e 5 ) 


2\/J + (x - X 0 ) + 2 (x-xo) + 

(15) 


The Hrst expansion is valid for regions to the left of x a , 
corresponding to the outer edge of the ESC and (at not 
too high a voltage) a quasi-LEN diffusion layer qualita¬ 
tively similar to Eqn. [3] The second expansion is valid 
within the ESC itself and can be integrated to find the 
voltage drop across the ESC to leading order, 

AVesc = -^E(-1-x 0 ) 3/2 (16) 

3e 

Note that while this gives a large (0(e -1 )) contribution 
to the voltage drop, hence dominating the response in 
terms of voltage, the lowest-order non-vanishing contri¬ 
bution to the concentration in the ESC arises from the 
first-order correction. It was shown in [19 ] that this ap¬ 
proach yields a well-defined minimum in the marginal 
stability of the ESC. Because this contribution affects the 
structure and resistance of the ESC and is also coupled 
to the selectivity of the membrane, we choose to employ 
it calculating the ESC-membrane Donnan condition. 

A rigorous treatment of the problem (see Ref. 0) 
involves defining appropriate boundary layer variables 
and careful consideration of the behavior of the (trans¬ 
formed) parameter x Q across the full range of applicable 
applied voltages in the limit e —> 0. Here we adopt a less 
rigorous approach, similar to that employed in 0 , to 


demonstrate the coupling between ESC structure, non¬ 
ideal selectivity, and the Ohmic-to-OLC transitional ESC 
resistance maximum 0. The idea is akin to solution¬ 
patching with a first-order correction and demanding 
that Eqns. [3] and [15] yield the same value for the con¬ 
centration at our estimated ESC edge, x = x 0 . 

Transforming to the boundary layer variables, E = 
(2 J + e) 1 / 3 GF and x — x 0 = (2 J + ) _1 / 3 e 2 / 3 z, the salt con¬ 
centration at the ESC edge is given by 

, ^ (e2J+) 2 / 3 G 2 F 2 (0) 

c(xo) = - - - (17) 

so that the ESC edge position can be approximated by, 
z 0 = ^-T-(2e) 2/3 G 2 ^^J“ 1/3 . (18) 


This requires the value of the E-field at x 0l or in terms 
of the boundary layer variables, the value of F(0), where 
F is the solution of the transformed Painleve equation, 

F" -^G 2 F 3 + zF = - 1 (19) 


with asymptotic solutions 



z < 0 
z > 0 


( 20 ) 


The solution is obtained numerically using a modified 
shooting method with the asymptotic solutions as bound¬ 
ary conditions. Initial guesses are supplied by higher- 
order corrections to the asymptotic expressions for the 
first few values of G. The solution for higher values is 
obtained by employing a (finite-difference) Taylor expan¬ 
sion of F(z; G) in SG. F{z ) is seen to decrease in the ESC 
with decreasing selectivity (see Fig. [2]. To use this result 
conveniently in subsequent calculations, the value of F(0) 
as a function of G is obtained using a fourth-order poly¬ 
nomial interpolation of the results. By evaluating the 
value of x 0 at the limiting salt flux density, it is shown in 
Fig. [3] that decreasing the selectivity increases the ESC 
thickness and that longer diffusion length to membrane 
thickness ratios are more sensitive to changes in the se¬ 
lectivity. A more pertinent question is how the selectivity 
changes couple to the current density itself. This requires 
use of Eqn. [9] to obtain the dependence of selectivity on 
applied current. A similar analysis was carried out in [Tlj 
up to the limiting current. Here we extend those consid¬ 
erations to the over-limiting case. 

As stated previously, we impose continuity of the 
electrochemical potentials, fi±(x) = lnc(a;) ± <fi(x) be¬ 
tween the ESC aX x = —1 and the membrane interior 
as a jump condition across the EDLs. This yields the 
following expression for the interior membrane boundary 
concentration at x = —1 in the under-limiting regime, 

CmM~ 1) = \/[1 - J + ( L - I)]" + N2 / 4 ( 21 ) 
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FIG. 2: Numerical solution to modified Paineleve equation 
lEan. flQll for several values of G (solid black). Large-z asymp¬ 
totes (Eqn. 1201) are shown also (dashed black.) The inset 
shows the interpolated value of E(0) as a function of G, used 
in subsequent calculations. 



FIG. 3: Estimate of the ESC width at the limiting salt flux 
density as a function of selectivity for e = 0.001. Two cases 
are shown, for diffusion layer length equal to the membrane 
thickness and for a diffusion layer much thicker than the mem¬ 
brane. For thinner diffusion layers, the result is (unsurpris¬ 
ingly) that the ESC is thinner overall and the change in thick¬ 
ness is less. 


contribution (such as the approach taken in M) , the 
qualitative picture would not change. 

The current-voltage relationship for the over-limiting 
case can be obtained simply by modifying the under¬ 
limiting (LEN) model fTTI.Il~2f . 


V = In 


Cm (1) ~|~-ZV/2 
Cm,«i(-l) + ^/2 


(G + 1) In 


CrG) 

c<(-1) 




(24) 


The above equation, like the internal membrane concen¬ 
trations in Eqns. [21] and [25] is obtained from the con¬ 
tinuity of electrochemical potentials. However, upon in¬ 
spection, it is clear that it may be thought of as a sum of 
voltage drops across the layers in the system. While it is 
generally true that one cannot arbitrarily add response 
coefficients in a multi-layer system [20], one can always 
add voltage drops themselves, provided such drops can 
actually be resolved in a straightforward way. Thus, one 
may insert the ESC voltage drop, Eqn. [TB] into the sum. 
It is worth emphasizing that Kirkwood’s argument re¬ 
garding the continuity of electrochemical potentials be¬ 
tween phases only applies between the membrane inte¬ 
rior and the interior of the ESC and across the EDL 
in the enriched side. Here, the EDLs are in fact thin, 
quasi-equilibrium boundary layers separating physically 
distinct phases. However this is not the case at the ar¬ 
tificial boundary, implied by x 0 . This is a mathemati¬ 
cal artifice corresponding to the poles in the asymptotic 
expansions, and not a Gibbs plane between physically 
distinct phases. So it is sufficient to demand continu¬ 
ity of concentration, as we did to obtain the estimate of 
x 0 . However, it is necessary to modify the concentration 
at the edge of the quasi-electro-neutral DL and the DL 
length accordingly. Doing so, one obtains the following 
approximate form for the current-voltage relationship in 
the over-limiting regime, 


V = In 


c m (l)+N/2 

)+N/2 


(G + 1) In 


Cr-(l) 

1 -J + (L+Xo) 


-CJ±Z(1) + AV ES c. (25) 


and in the over-limiting regime, using the concentration 
c(x 0 ) of Eqn. [T71 


Cm,ol( 1) 



( 22 ) 


and the usual expression for the concentration at the en¬ 
riched side, 

c m ( 1) = \]{J+(L-l) + lf + N*/A (23) 

To be consistent with the estimate of the ESC width, we 
have kept the first-order correction term, even though 
this amounts to a change in the membrane interior con¬ 
centration of order 0(e773)_ Were we to neglect this 


III. RESULTS 

The dependence of the selectivity on salt flux density 
may be extracted directly from Eqn. [9] by substituting 
the concentration at x = —1, using Eqn. [Til or [551 and 
that at x = 1 Eqn. 1231 for the appropriate regime. This 
is shown in Fig. [I] The ratio of the absolute value of co- 
to counter-ion flux density is plotted against overall salt 
flux density for a range of parameters, with an increase in 
the ratio corresponding to loss of selectivity. The range 
of salt flux density is the same for all three cases, from 
zero to five times the limiting value. The selectivity de¬ 
creases as the salt-flux is increased because the difference 
in co- and counter-ion electrochemical potential across 
the charge-selective region increases. This corresponds 
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FIG. 4: Selectivity variation with salt flux density in terms 
of ratio of co- to counter-ion flux density for various system 
parameters where the standard case is IV = 10, L = 3, and 
T>m = 1. Each plot shows a comparison of the standard to 
other cases; (a) varying fixed charge density N, (b) varying 
diffusion length L, and (c) varying membrane diffusivity, D m . 


to both an increase in conductivity and increasing co¬ 
ion concentration within the charge-selective region. In 
agreement with previous studies 11], the change in the 
under-limiting regime is rather modest. However, since 
the enriched concentration has no imposed upper bound, 
this difference can be expected to increase as J+ exceeds 
■ 

In Fig. 0Jr, the change in selectivity is shown for 


several values of the fixed charge density. Unsurpris¬ 
ingly, as the charge density increases-corresponding to 
increased restriction of the co-ion flux-both the over¬ 
all selectivity and change in selectivity decrease. The 
cases in Figs QJ) and c bear some relationship to the ef¬ 
fects of field-focusing in the sense that increased field- 
focusing corresponds to both increasing flux density in 
the focused region and effective shortening of the diffu¬ 
sion length [21]. For different diffusion lengths L , the 
total change in selectivity across the range of salt flux 
density, from J + = 0 to J+ = 5 J+,u m , is the same for all 
L. However, the overall selectivity itself decreases with 
increasing L because the concentration at the enriched 
boundary is higher for the same salt flux density. The 
same is essentially true of varying T> m ] the total change 
in selectivity is independent of the effective diffusivity in 
the charge-selective region. However, while the limiting 
salt flux density in the diffusion layers is the same for all 
cases, a decrease in effective diffusivity may be thought 
of as an increase in effective salt flux density inside the 
charge-selective region. So a comparably small change 
in concentration at the enriched side, corresponding to a 
comparably small total change in electrochemical poten¬ 
tial across the membrane, leads to lower overall selectiv¬ 
ity without shifting the voltage. 

To validate our approach, we compare the results of 
the above model with numerical calculations made us¬ 
ing COMSOL software, shown in Fig. [5] The agreement 
between the approximate analytic model and the numer¬ 
ical model is quite good, diverging appreciably only at 
the highest applied voltages where the analytic model 
under-predicts the current density by a few percent. The 
two models yield essentially identical selectivity response 
across the entire voltage range, showing that the bulk of 
the selectivity change occurs at higher voltages. 

Concerning the resistance maximum, both the maxi¬ 
mum value of the resistance and the voltage at which it 
occurs are related to the changes in selectivity, and thus 
tied to the other parameters in the system. The change 
in selectivity itself does not require an accurate value of 
the concentration in the ESC; similar change would be 
obtained by simply equating the over-limiting concentra¬ 
tion to zero. However, the very existence of the maxi¬ 
mum itself arises from the fact that not only is the con¬ 
centration non-zero, but increases with current beyond 
the limiting value (see Eqn. [13. The basic mechanism 
behind the maximum may be understood as follows: Ini¬ 
tially, the increase in the extent of the depleted region, 
as indicated by x a , causes a rapid increase in the ESC 
resistance. However, as the ESC concentration increases, 
the effective conductivity of the depleted region increases 
both within the ESC and at the edge of the quasi-electro- 
neutral region, giving rise to the maximum. 

By considering the variation of the concentration 
c( x 0 ) (Eqn. 1171) along with selectivity variation one can 
understand the results shown in Fig. [6] The parameters 
correspond to the same cases as shown in Fig. [3 Both 
total resistance and ESC resistance are shown for each 
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FIG. 5: Comparison of theoretical model (dashed black) to 
COMSOL calculations (solid blue) for a ID PNP model, both 
using the standard parameter set N = 10 ,L = 3, and D m = 1. 
(a) Current-voltage response is shown in the top figure and 
changes in selectivity in the bottom (b). Inset in top figure 
shows discrepancy between approximate model and COMSOL 
grows with voltage. 


case. For the range of parameters considered here, the 
resistance is almost totally dominated by the ESC resis¬ 
tance. 

Fig. [S}i shows the resistance-voltage curves for differ¬ 
ent values of the fixed charge density, N. For large val¬ 
ues of IV, the selectivity is high and, as seen in Fig. @Ji, 
changes little with applied voltage. The co-ion concen¬ 
tration inside the charge-selective region remains low and 
the counter-ion concentration is governed primarily by 
the value of N itself. Thus G remains very close to the 
ideal value, G = 1 for all voltages and the ESC con¬ 
centration increases slowly with voltage, being governed 
primarily by the slow increase of J/ . So the overall 
resistance is high and the slow change in concentration 
imply a maximum resistance at higher voltage than the 
other cases. For smaller values of IV, the change in G 
(equivalently, the ratio |j_/j + |) is larger, resulting in in¬ 
creasing G faster than the decrease in F( 0). Hence there 
is a larger increase in the ESC concentration and lower 
resistance occurring at a lower voltage for smaller N. 

The dependence of the resistance-voltage curve on dif- 


(a) 



(c) 



FIG. 6: Resistance vs. voltage curves, for cases matching 
figure 4. The ‘standard’ case (solid red) has N = 10, L — 
3, and D m . The total resistance (dashed black) and ESC 
resistance (solid blue or red) are shown for each case. 


fusion length L is more indicative of the changes in the 
limiting salt flux density and in the enriched concentra¬ 
tion than the increased resistance of the longer diffusion 
layers. This can be expected to differ from a true 2D/3D 
case, where field-focusing plays a strong role [2l| , or a sit¬ 
uation in which non-uniform fluid flow can alter the de¬ 
pleted concentration profile ((_>. 7l. The maximum occurs 
at higher voltages for longer diffusion lengths because the 
total change in selectivity occurs over a smaller range of 
salt flux density while the actual value of the salt flux 
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FIG. 7: COMSOL geometry sketch, note that L here does 
not correspond directly to that in the ID model. 


density is lower for longer L (see the discussion associ¬ 
ated with Fig. 0, corresponding to a smaller increase in 
ESC concentration for a given change in total applied 
voltage. Thus it takes a larger overall voltage to reach 
the resistance maximum. The reasoning for changing dif- 
fusivity is somewhat similar in the sense that higher ef¬ 
fective salt flux in the membrane for smaller values of T> m 
corresponds to an overall higher selectivity, hence lower 
co-ion concentration and flux. However, because it does 
not correspond to an actual modification of the concen¬ 
trations or electric potentials adjacent to the membrane, 
the voltage does not shift. 

To gain some insight into the influence of field- 
focusing on the structure of the ESC, selectivity varia¬ 
tion, and dc response, a 2D numerical solution for the 
PNP equations is obtained using COMSOL. The geome¬ 
try is shown in Fig. [7] A uniform volumetric charge den¬ 
sity, N = 10, is specified in the constricted ‘nanochan- 
nel’ region. No surface charge is specified on the channel 
walls, so that we may study the effects of field-focusing 
on the ESC in isolation from surface-conduction. Fig.[51i 
shows the current-voltage characteristics in the main fig¬ 
ure, which are qualitatively similar to those obtained for 
either decreasing diffusion-length L or decreasing T> m in 
the ID model. However, note the selectivity response 
(inset) shows that the increased field-focusing results in 
both higher selectivity initially, and a greater overall 
change in selectivity with increasing salt flux (equiva¬ 
lently, voltage.) It is important to keep in mind that the 
numerical calculations are done under potentiostatic con¬ 
ditions, as opposed to the galvanostatic conditions em¬ 
ployed in analytic model. Therefore, the voltage range 
is kept fixed and the salt-flux and current densities are 
determined in response to the voltage. Thus the same 
voltage range produces a change in selectivity over a 
smaller range of salt-flux density for an increasing de¬ 
gree of field-focusing. The resistance-voltage curve in 
Fig. [8 Jd further illustrates the difference. While the value 
of the resistance increases, in this case corresponding to 
a smaller nanochannel height for a fixed microchannel 
height, the resistance maximum shifts to a higher volt¬ 
age. This stands in contrast to both decreasing V m and 
decreasing L, where by analogy we would expect the in¬ 
creased degree of field focusing to correspond to either no 
shift in the voltage of the maximum, or a shift to lower 


(a) 



(b) 



FIG. 8: Comsol cases L = 2 ,N = 10, bulk diffusivity ev¬ 
erywhere and (i) H = h = 1 corresponding to the standard 
case of the analytic model, (ii) H = 1 ,h = 1/2, (iii) H = 1, 
h = 1/4. (a) I-V with selectivity ratio inset, (b) R-V 

voltages. 

To better understand what is happening in the 2D 
case, we consider the evolution of the 2D ESC with 
voltage. The centerline salt concentration near the 
nanochannel entrance is plotted in Fig. [9] for voltages 
below and above the maximum. Note that, as one might 
expect, increasing field-focusing corresponds to a smaller 
ESC at the same voltage. However, the ESC concentra¬ 
tion undergoes less of an increase with voltage when the 
degree of field-focusing is higher. While this seems at first 
counter-intuitive, it is worth noting that the range in salt 
flux density for the same voltage range is reduced consid¬ 
erably. Thus, at a given voltage, a lower salt flux density 
drives the change in ESC salt concentration leading to a 
slower increase in ESC concentration with increased field 
focusing. 


IV. CONCLUSIONS 

We have presented theoretical and numerical results 
concerning the connection between variation of selectiv¬ 
ity, ESC structure, and dc over-limiting resistance of a 
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FIG. 9: Depleted salt concentration profile for numerical 2D 
case, evaluated on centerline near nanochannel entrance with 
L = 2,N = 10, bulk diffusivity everywhere. Solid lines 
evaluated for V = 1000 and dashed lines for V = 7000 (i) 
H = h = 1 corresponding to the standard case of the analytic 
model, (ii) H = 1 ,h = 1/2, (iii) H = 1, h = 1/4. 

3-layer model electro-diffusive system with a non-ideal 
charge-selective element. This approach extends previ¬ 
ous studies concerning the change in selectivity for a sim¬ 
ilar system in the under-limiting regime 0 - More im¬ 
portantly, the present study represents progress towards 


developing a model capable of capturing the OLC dc re¬ 
sponse of a fabricated micro-nanochannel system. The 
majority of previous work has relied on models which as¬ 
sumed ideal permselective transport, which make it dif¬ 
ficult to understand the relationship and interaction of 
the ESC with other OLC mechanisms such as surface 
conduction fU, ,22j or non-uniform fluid flow @, 0. It 
was recently shown that non-ideal selectivity can lead 
to a bulk electro-convective instability within the locally 
electro-neutral approximation [23] , and the present work 
serves to demonstrate that the structure of the ESC itself 
is affected by both non-ideality and field-focusing. It is 
reasonable to assume such structural change can affect 
electro-convective stability as well, and inclusion of these 
effects could be expected to enrich the present under¬ 
standing of electro-convection in OLC and high-current 
transport in experimental micro-nanochannel devices. 
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